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ABSTRACT 

It is shown that the matrix models which give non-perturbative definitions of string and 
M theory may be interpreted as non-local hidden variables theories in which the quantum 
observables are the eigenvalues of the matrices while their entries are the non-local hidden 
variables. This is shown by studying the bosonic matrix model at finite temperature, with 
T taken to scale as 1/N. For large N the eigenvalues of the matrices undergo Brownian 
motion due to the interaction of the diagonal elements with the off diagonal elements, giving 
rise to a diffusion constant that remains finite as — > oo. The resulting probability density 
and current for the eigenvalues are then found to evolve in agreement with the Schroedinger 
equation, to leading order in 1/N, with h proportional to the thermal diffusion constant 
for the eigenvalues. The quantum fluctuations and uncertainties in the eigenvalues are then 
consequences of ordinary statistical fluctuations in the values of the off-diagonal matrix 
elements. Furthermore, this formulation of the quantum theory is background independent, 
as the definition of the thermal ensemble makes no use of a particular classical solution. The 
derivation relies on Nelson's stochastic formulation of quantum theory, which is expressed in 
terms of a variational principle. 
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1 Introduction 



In this paper we describe a new proposal concerning the relationship between general relativ- 
ity and the quantum theory. This is that matrix formulations of string, or M. theory [Q, 0, 
|J, which are known to reproduce general relativity, at least perturbatively, are also hidden 
variables theories. More precisely, given the dynamics postulated by these models, we show 
that the finite temperature classical statistical mechanics of the matrix elements reproduces 
the quantum theory of the matrix eigenvalues. This happens when the temperature is scaled 
appropriately, as T pa 1/N as N is taken to infinity. Thus, these theories contain not only a 
non-perturbative definition of string theory, they contain a reformulation of quantum theory 
in terms of the ordinary statistical mechanics of a set of non-local variables which are the 
matrix elements. 

By a hidden variables theory we mean generally a theory which purports to give a more 
detailed description of individual events and processes for which the quantum mechanics 
yields only probabilistic predictions. It is thus a theory whose ordinary statistical mechanics 
reproduces, at least to a certain order of approximation, the conventional quantum me- 
chanics. The name hidden variables refers to the existence of degrees of freedom beyond 
the observables of the quantum theory, whose statistical fluctuations are the source of the 
quantum fluctuations and uncertainties of quantum theory. 

We know from the observed violations of the Bell inequalities that any viable hidden 
variable theory should be non-local. This leads to a natural suggestion, which is that the 
non-local hidden variables describe relationships between the local degrees of freedom ||. A 
simple hypothesis is that there is such a relational hidden variable for each pair of degrees of 
freedom of a quantum theory. This then suggests formulating a non-local hidden variables 
theory in which the fundamental, hidden degrees of freedom are elements of matrices while 
the observables are the eigenvalues of the matrices. 

Such a hidden variables theory was developed some years ago and shown to reproduce 
the predictions of non-relativistic quantum mechanics for an iV-body system]^]. When the 
elements of an N x N matrix are put into a thermal bath, it was found that under certain 
conditions, the probability density and current for the eigenvalues behaved to leading order 
in 1/yN, in a way that reproduced the iV-body Schroedinger equation. Other non-local 
hidden variables have been formulated and studied, for example, the Bohm model ||. 

Thus, it is not very difficult to make non-local hidden variables theories. What is more 
challenging is to invent such a theory that solves other problems besides the issues in the 
foundations of quantum mechanics. Since a hidden variables theory must be non-local, it 
seems very likely that a true one would have something to do with the structure of space 
and time, and hence with quantum gravity^} 

It is then remarkable that matrix models which make use of essentially the same idea 
as the hidden variables theory described in||-that the real physical degrees of freedom are 
matrix elements, while the eigenvalues correspond to the observables-were developed as rep- 

*A different proposal to unify relativity with quantum theory in the context of a hidden variables theory 
has been proposed by 't Hooft||. 
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resentations of string and M. theory. These matrix formulations of string theory have solved 
certain problems, but during their study a certain puzzle has emerged concerning their quan- 
tization. 

This is a consequence of the fact that one defines string theories from compactifications 
of the matrix models. To define a quantum string theory one first picks a classical solution of 
the matrix model, which gives information about the background fields including the space 
time geometry. One then proceeds to define a quantum theory with respect to that classical 
background. 

Of course, one can define formal quantizations, for example by means of path integral, 
which appear to be background independent. However as soon as one attempts to define the 
quantum theory precisely, for example, in terms of a BRST quantization, one finds that the 
precise definitions of the Hilbert space and BRST charge Q do depend on the background. 
Or, if one wants to define the quantum theory precisely in a path integral context, one faces 
the problem that the Lorentzian path integral does not appear to be well defined. However, 
the usual solution to this problem, which is to make a Euclidean continuation, is itself back- 
ground dependent, because it depends on a choice of time in a particular classical spacetime 
manifold. No background independent definition of continuation between Euclidean and 
Lorentzian theories is known. 

This problem, of the background dependence of the quantizations, is a major problem 
for string and Ai theory. As the theory classically has many solutions, which define different 
spacetime backgrounds, it seems that there ought to be a background independent way to de- 
fine the quantum theory. If the fundamental theory is a quantum theory it should be that the 
different quantizations defined around the different backgrounds should be approximations 
to a single exact quantum theory. 

This is part of the motivation for the present proposal. Perhaps the background inde- 
pendent theory is not a conventional quantum theory, but some deeper theory, which can 
be approximated by conventional quantum theories when the state defines a fixed spacetime 
background. Such a theory might be approached in different ways, but here we investigate 
the hypothesis that it might be a hidden variables theory. As we will show, the hint from 
earlier work is correct, and a matrix model of the kind studied in string and M. theory can 
serve as a non-local hidden variables theory which can then reproduce a quantum theory for 
the eigenvalues of the matrices to leading order in 1/N. 

That is, rather than quantizing the classical matrix model in some conventional fashion, 
we will simply assume that the off diagonal elements of the matrix model are in a classical 
thermal state. We will find that the quantum theory for the eigenvalues can be reproduced so 
long as the temperature is scaled in a certain way with N. This formulation of the quantum 
theory is by definition background dependent, because the definition of the thermal ensemble 
makes no reference to any particular classical solution. 

To see how this happens it is useful to consider the diagonal elements of the matrices, 
which become increasingly free at low temperatures, as analogous to the pollen grains in 
classical Brownian motion. The off diagonal elements are then analogous to the molecules 
whose constant collisions with the grains cause them to move with a Brownian motion. 
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Indeed, the diagonal elements are subject to random forces from their interactions with the 
off-diagonal elements. The off diagonal elements are small at low temperature, but as we 
increase N their effects on the diagonal elements are greater due to their greater number. 
The result is that the interaction of a diagonal elements with a large number of off-diagonal 
elements introduces a Brownian motion, which is transferred to a Brownian motion of the 
eigenvalues at low temperature. That is, the randomness of the local variables-the diagonal 
elements and the eigenvalues is due to their interactions with a much larger number of 
non-local variables. 

To find interesting behavior we have to scale T in an appropriate way with iV as we 
take the former to zero and the latter to infinity. In fact, we find that the model behaves 
critically when we scale the temperature so that T ~ 1/N. In this case the off diagonal matrix 
elements are of order 1/yN. However their collective effects on the diagonal elements, and 
hence the eigenvalues, remain as N — > oo. One such effect is that the diffusion constant 
which measures the Brownian motion of the eigenvalues remains finite as N — > oo. 

Of course, the idea that quantum statistics might just be ordinary statistics in an unusual 
context is an old one. In particular, Nelson|7| has proposed a stochastic formulation of 
quantum theory according to which the quantum description of a particle is derived by 
modifying the classical description solely by the addition of a universal Brownian motion, 
which satisfies certain special properties. Chief amongst them is that the Brownian motion 
is non-dissipative, in that energy and momentum are still conserved. Nelson shows in |7j] that 
when a classical particle is subject to such a non-dissipative Brownian motion, its probability 
density and current evolve in a way which is equivalent to that given by the Schroedinger 
equation. 

Nelson's formulation plays a key role in the present work, in that we show that the 
stochastic formulation of quantum theory is recovered to leading order in 1/N, for the eigen- 
values of the matrices in our model[|. 

In the present paper we study a bosonic matrix model, and show that it is indeed a non- 
local hidden variables theory. The extensions of this work to the supersymmetric matrix 



models associated with string and M. theory are in progress with Stephon Alexander |L0 . 

In the next section we describe the matrix model we will study. A variational principle 
related to Nelson's formulation of quantum theory is presented in section 3, where the basics 
of the theory of Brownian motion are reviewed for those unfamiliar with it. In section 4 
we describe the basic physical picture which suggests the connection between the classical 
statistical mechanics of the matrix model and quantum theory. In section 5 we estimate the 
dependence of the relevant diffusion constants on N , T and other parameters. Finally, the 
derivation of the Schroedinger equation for the eigenvalues is given in section 6. 



2 Note that one can quantize the free bosonic string directly using Nelson's stochastic quantum theory |ll[ 



5 



2 The model 



We study a bosonic matrix model which is the bosonic part of the models used in string and 
M. theory [fll [3], M. The degrees of freedom are d N x N real symmetric matrices Xj, with 
a = 1, d and i, j = 1, N. The action is, 

S = n | rftTr [x 2 + a; 2 [X a , X 6 ] [X a , X fe 

We choose the matrices X a to be dimensionless. u is a frequency and \i has dimensions of 
mass • length 2 . We do not assume % = 1, in fact, as we aim to derive quantum mechanics 
from a more fundamental theory, fi is not yet meaningful. We will introduce H as a function 
of the parameters of the theory when we derive the Schroedinger equation as an approximate 
evolution law. We may note that the parameters of the theory define an energy e = [iio 2 . 
The basic idea is that the off diagonal matrix elements of X a will be the non-local hidden 
variables. The physical observables will be defined to be the eigenvalues A" of the matrices. 
We will put the system at a small, but finite temperature, the result of which will be that 
the matrix elements undergo Brownian motion as they oscillate in the potential. It follows 
from linear algebra that the eigenvalues also undergo Brownian motion. We will see that 
the parameters of the theory can be scaled with N in such a way that Nelson's stochastic 
formulation of quantum mechanics is realized for the eigenvalues. When T = the matrices 
must commute with each other so as to achieve the vanishing of the potential energy. This 
means that it is possible to simultaneously diagonalize them. When T is finite, but small 
compared to e, the off diagonal elements will on average be small. As a result, it is useful to 
split the matrices into diagonal and off-diagonal pieces, 



(2) 



where D a = diagonal^", d%) is diagonal and Qj- has no diagonal elements. Since the Qj 
are dimensionless we will expect them to scale like a power of T/fiuj 2 . We then write the 
action [1] as 



S 



The theory of the cf s alone is free, 



dt 



c d + c Q + c 



int 



at 

while the theory of the Q's alone has the same quartic interaction 



/i 



J2m) 2 +co 2 [Q a ,Q b ][Q a ,Q b 



aij 



(3) 
(4) 

(5) 



The interaction terms between the diagonal and off-diagonal elements are 
C mt = 2^ 2 [-« - d«f{Q£f - ^ - d«){d\ - 4)QjQ 6 } + 2K - d?)<2?[Q a , Q 

abij 



(6) 



6 



We note that the model has a translation symmetry given by 



d? — d? + t/ 1 . 



(7) 



The result is that the center of mass momentum of the system is conserved 

3 The statistical variational principle 

We now introduce the basic ideas behind the stochastic formulation of quantum mechanics 
of Edward Nelson. We introduce Nelson's formulation by means of a variational principle, 
which we call the statistical variational principle. This is closely related to variational prin- 
ciples previously introduced by Guerra||| and by Nelson 0. We then will formulate it in 
hamiltonian langauge, which will provide insight into how the linearity of the state space of 
quantum theory emerges from the theory of Brownian motion. 

3.1 S ensembles 

Consider a dynamical system living on an n-dimensional configuration space coordinatized 
by x a . The dynamics can be described in terms of a Hamilton- Jacobi function, S(x, t) which 
solves the Hamilton- Jacobi equation 



A particular solution S defines a family of classical trajectories whose momenta at any point 
x a are defined by 



There are many solutions to the Hamilton- Jacobi equation, each of which defines a congru- 
ence of classical trajectories. A statistical description of the system may be given in terms 
of a probability density p(x,t) and a probability current v a (x, t), which together satisfy 



Since the probability is conserved we may always assume / d n xp(x,t) = 1. 

Now we will do something unusual. Let us restrict attention to an ensemble of classical 
trajectories whose evolution is determined by a particular solution S of the Hamilton- Jacobi 
equation. We may call this an S-ensemble. These have a probability density p$(x,t). Since 
the momentum is determined by S, so must be the probability current. We then have, 



s+^-(d a sr+u=o 





(9) 




(10) 




(11) 



so that the probability conservation equation is now 



p+-d a (pd a S) = 



(12) 
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The restriction of attention to a statistical ensemble with fixed solution S is unusual, but it 
does not take us out of the domain of classical physics. The total probability density may 
be recovered formally as 

Ptotal=^2PS (13) 

s 

where the sum is over all solutions to the Hamilton- Jacobi function. Of course, the probabil- 
ity current does not add. Now we may notice that such ensembles are given by the solutions 
of a simple variational principle 



I[p, S} = J dt J d n xp(x,t) 



S + ^-(d a S) 2 + U(x) 
2m 



(14) 



The equations that arise from varying p and S are, respectively, the Hamilton- Jacobi equation 
(§) and the probability conservation equation (Jl2|). Thus, the variational principle describes 
an ensemble of trajectories, each of which evolves according to the Hamilton- Jacobi equation, 
so that the current velocity is proportional to d a S 

We note that the action and equations of motion are invariant under time reversal with 

t -> -t, S^-S, p^ p (15) 



3.2 A brief review of the theory of Brownian motion 

Nelson's stochastic version of quantum theory may be formulated in the language we have 
just introduced. To do this we assume that in addition to the classical motion, the particles in 
our ensemble are subject to a Brownian motion. This Brownian motion is, however, unusual, 
in that it does not alter the condition that each trajectory in the ensemble is governed by 
the same solution of the Hamilton- Jacobi equation. We will see that this requirement can be 
met, by altering in a small way the Hamilton- Jacobi equation itself, to take into account the 
fact that trajectories undergoing Brownian motion are not differentiable. Nelson calls the 
resulting Brownian motion dissipationless Brownian motion, as energy and momentum are 
still conserved. To describe this dissipationless Brownian motion we may use the language 
of stochastic differential equations^ In this language the small change in time of a trajectory 
is given by 

Dx a = b a (x(t),t)dt + Aw a dt>0 (16) 
for small positive changes in time and 

D*x a = -b* a {x(t),t)dt + Aw* a dt<0 (17) 

for small negative changes in time. b a and b* a are called the forward and backwards drift 
velocities. They describe the average motion of the particles in the ensemble. They are 
defined by 

,„, , , x a (t + At) - x a (t) . , 

b a (x, t) = hm < ^ ±± > x(t)=x (18) 



3 For reviews see 
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and 

,*„/ s , x a (t) - x a (t - At) , . 

b* a (x, t) = Hm < u ^ l i > x{t)=x (19) 

The different elements of the ensemble are distinguished by their Brownian motion, which 
is given by a Markov process defined by 

< Aw a Aw b >=- < Aw* a Aw* b >= vdtq ab (20) 

and 

< Aw a Aw* b >= (21) 

Here q ab is a metric on the configuration space and v is the diffusion constant. The averages 
< ... > are defined with respect to the ensemble. Thus, for any function F(x) on the phase 
space 

< F >= J d n xp(x,t)F(x) (22) 

From these basic definitions one can derive easily the forward and backwards Fokker-Planck 
equations 

p = -d a {pb a ) + vV 2 p (23) 
p = -d a (pb* a ) - uV 2 p (24) 
From these the current conservation equation follows with 

v a = ^{b a + b* a ) (25) 

The difference between the forward and backward drift velocities is called the osmotic ve- 
locity. From the Fokker-Planck equation it satisfies 

u a = i(6 a -6* a ) = vd a \xip (26) 

We thus see that the diffusion constant measures the extent to which the paths are non- 
differentiable, so that the forward and backwards drift velocities are not equal. This is of 



course possible because they are defined in eqs. ( |18[|19D in such a way that the limit At — > 
is taken after averaging over the ensemble. 

3.3 Quantum Brownian motion and Nelson's stochastic formula- 
tion of quantum theory 

With this quick survey of Brownian motion over, we return to the case of interest, which is 
an ensemble of trajectories that share the same Hamilton- Jacobi function, S. We want to 



preserve the property (11) that the current velocity is proportional to the gradient of the 
Hamilton- Jacobi function, but we want to find a way to include Brownian motion within the 
ensemble defined by a particular Hamilton- Jacobi function. One way to approach this is to 
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modify the statistical variational principle to include the effects of Brownian motion. It is 
not hard to see that this is possible, and that the right extension of the variational principle 
is 



I v \p,S\ = j dt J d n xp(x,t) 



1 771 7/^ 

S + ^(d a S) 2 + — (d a hip) 2 + U(x) 



(27) 



To see why let us use the fact that d a S is proportional to the momentum. Thus we have for 
smooth motion, on solutions to the variational principle, 



[dt [ d n xp(x,t) — (d a S) 2 = I dt [ d n xp(x,t) — (p a ) 2 
J J 2m J J 2m 

m 



dt J d n xp(x,t)j(x a ) 2 

, f m / ^ m i f X a (t + At) - X a (t)Y . , 

it j cfVfz, ()- km (28) 

However, for a Brownian motion process, the limit in the last line is not defined. So this is 
not a consistent variational principle when v ^ 0. To define a variational principle that is 
well defined for the case of Brownian motion where the paths are non-differentiable we need 
to take the limit that defines the derivative outside of the ensemble average. Thus, we may 
define instead, 

Jdtj d«xp(x,t)±( P r = J dthmj d n xp(x,t)^ Mt + ^-Afl V (29) 

This form of the integrand appears, however, to lack invariance under time reversals, eq. 
(15). This is because when the paths are non-differentiable b a and b* a may not be equal. 
However, we may notice that under the time integral we can write 

m f x a (t + At) - x a (t)\ 2 r r m f x a (t + At) - x a (tY 



. , f x a (t) - x a (t - At)\ 
+P (x,t-At) l ) 



At 

2 



(30) 



Now, we notice that 

p(x, t - At) = p(x, t) - Atp(x, t) (31) 

As p(x, t) is given by the Fokker Planck equation the second term leads to terms that are 
well behaved and vanish as At — * 0. Thus, we can take the average and then the limit, using 
(0,0) to find, 

/ dthraj d"xp(x,t)^ ^(t + ^-^ J = J dt J d "xp(x,t)^m 2 + (bn 2 + C] 

dt J d n xp(x,t)^[(v a ) 2 + (u a ) 2 + C] 

(32) 
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Here C is an infinite constant, which is equal to 

C = ud lim — (33) 

Thus, we have found that we can extract an infinite constant from the action, leaving us 
with a finite piece that is well defined and time reversal invariant. So we have, 

/* & /^M)f plMzgW )' = /,ft/ ( Pxp(x, t )[^(4S)» + ^(«.lnp)»l + ^ 

(34) 

Constants, even infinite constants, play no role in classical action principles. Hence, putting 
this last result back in the definition of the action principle we see that the effect of exchanging 
the order of the integral and the limit is the new action principle (0). Thus we see that while 
we can have a Brownian motion within the trajectories defined by a given Hamilton- Jacobi 
functional, the definition of Brownian motion requires that we modify the Hamilton- Jacobi 
equation, so that we have a consistent variational principle, even in the presence of non- 
differentiable paths. 

From the argument just given we see that the variational principle (|27|) is equivalent to 
defining the time derivatives so that they are taken outside the ensemble averages. This 
yields the standard variational principle when the trajectories are smooth, because then 
it doesn't matter in which order we take the limits involved in the definition of the time 
derivatives and the ensemble average. But for non-differentiable paths the order matters 
and we see that we must take the limit defining the time derivative after that defining the 
ensemble average. As we have seen, up to an infinite constant which may be ignored, this is 
equivalent to adding the new term ^^^-(daln p) 2 to the variational principle. 

The new Hamilton- Jacobi equation follows from (|27|) by varying by p. We find that 

S+^-(d a S) 2 -2mu 2 ^-\/ 2 ^+U = (35) 
2m y/p 

The Hamilton- Jacobi functional is thus modified by a new potential term which is a function 
of the probability density. We recall that this unusual feature follows because it is the 
gradient of the probability density that measures the importance of the non-differentiability 
of the paths. The probability conservation, however is not modified. Now, the big surprise 
is that (|35| ) and (|T2| ) are nothing else than the real and imaginary parts of the Schroedinger 
equation, with 

V{x,t) = ^Jp{x,t)e lS/h (36) 

and with 

H = 2vm. (37) 
Thus, the variational principle (|27j), in the presence of Brownian motion is equivalent to 



ih 



dt 



2 



tf(z,t) (38) 
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One way to understand this is the following. From the point of view presented here, a 
quantum state is nothing more nor less than an ordinary statistical ensemble of Brownian 
motion trajectories, which share a single hamilton-jacobi function, S, where that S is itself 
a solution to a modified Hamilton- Jacobi equation, modified to take into account the change 
in the definition of the momentum necessary when the motion is Brownian. This is the basic 
message of Nelson's stochastic quantum theory. 



3.4 Hamiltonian formulation of the statistical variational principle 



can 



To find the hamiltonian formulation of the statistical variational principle note that 
be written 

r [p, S}= J dt J d n x [Sp - H[p, S, x}} (39) 
where the Hamiltonian density is 



H[p, S,x]= p 



I 777 //^ 

— (d a S) 2 + —(d a \np) 2 + U(x) 



(40) 



Thus we see that the probability density p can be considered as a conjugate coordinate with 
S its conjugate momentum, so that we have an infinite dimensional phase space with 



{p(x),S(x')} = 5 n (x',x) 



The Hamiltonian, 



H 



d n xH 



(41) 



(42) 



is then conserved in time. It is easy to check that Hamilton's equations of motion are the 
Hamilton- Jacobi equation, ( |55D and the probability conservation equation fll2"|) . We note 
that this is true for any value of the diffusion constant v so that this is true in both classical 
and quantum theory. To get more insight into how the linearity of quantum theory has 
emerged from the theory of Brownian motion, we can write out the conserved hamiltonian, 



H 



d n xp 



I 777 //^ 

— (d a S) 2 + —(d a \np) 2 + U(x) 



(43) 



This seems a very non-linear expression, but it is nothing but the expectation value of a 
linear operator. To see this we rewrite it slightly as 



-iS/h 



2m 



H = J d n Xyfpe 

Using (|36| ) and fl3~7|) this is easily seen to be equal to 

H = f d n xm^ 



VP 



+ U{x) 



iS/h 



(44) 



(45) 
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with 

H = -^-V 2 + U(x). (46) 

Thus, the very non-linear seeming equation H = is seen to be actually equivalent to the 
linear Schroedinger equation. We further see that the conserved Hamiltonian which arises 
from the statistical variational principle is exactly equal to the expectation value of the 
Hamiltonian operator in the quantum theory. Thus, the conservation of the Hamiltonian in 
the statistical variational principle is equivalent to the conservation of the expectation value 
of the Hamiltonian operator in the quantum theory. 



4 The physical picture 

We will now return to the matrix model we introduced in section 2. We will show that 
the ordinary statistical physics of this model has a critical behavior when the off diagonal 
sector is heated to finite temperature and the large N limit is taken with the temperature 
scaled so that T w 1/N. We will further see that a feature of the critical behavior in this 
limit is to reproduce quantum mechanics. That is, to leading order in 1/N the evolution 
of the probability density and current for the eigenvalues of the matrixes is equivalent to 
that given by the free Schroedinger equation. To show this we will apply what we have 
learned in the last section to the matrix model. We will formulation an S ensemble for 
the matrix model in terms of the statistical variational principle. At the fundamental level 
the dynamics is formulated in terms of the variational principle without Brownian motion, 
i.e. the simple variational principle (|I4]) . Thus, the whole matrix model is in an ordinary 
statistical ensemble. In the next section we will study the behavior the matrix model at low 
temperature and large N. We see that when we pick T « 1/N the off diagonal elements 
scale as ly/N ; this makes sense as they must go to zero at T = 0. We also see that to leading 
order the off diagonal elements can be seen to move harmonically in an average field given 
by the average values of all the other off-diagonal elements. The diagonal elements are not 
required to vanish as T — > 0, so they remain of order unity. However, the diagonal elements 
move in a random potential given by the oscillations of all the off diagonal elements. The 
result is that the diagonal elements pick up a random Brownian motion, on top of their free 
motion. This Brownian motion is then also experienced by the eigenvalues. We will see that 
when the model is scaled critically, the diffusion constants for the diagonal elements and 
eigenvalues go to constant limits as iV — > oo and T — > 0. We then want to study the effect of 
the Brownian motion of the eigenvalues. To do this we derive an effective statistical action 
for the probability distribution of the eigenvalues by averaging the statistical variational 
principle for the whole model over the values of the matrix elements. This is the task of 
section 6. We see that the Brownian motion term in (P7j ) emerges naturally as a term 
in the effective statistical action for the eigenvalues, as a result of the induced Brownian 
motion just described. Furthermore, in that limit the conserved energy of the variational 
principle of the whole model reduces to the conservation of the expectation value of the free 
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hamiltonian operator for the eigenvalues. Thus, in the large N limit Nelson's stochastic 
formulation of quantum theory emerges naturally as a description of the statistical behavior 
of the eigenvalues. 



5 Estimates for the diffusion constants at low temper- 
ature 

In this section we thus investigate the consequences of putting the matrix model in a thermal 
bath at a temperature T. We have two tasks. The first is to understand how various 
quantities of interest scale with T and N, at low temperatures. The second is to derive 
estimates for the diffusion constants for the matrix elements and eigenvalues that are good 
for low temperature and large N. By low temperature, T, we will mean that the ratio T / ^iuo 2 
is small. It will be convenient to scale this ratio with N, so we define 

T t 

(47) 



8(rf-l)/icj 2 Np 

with t dimensionless and p a power. The factor of 8(d — 1) is inserted for later convenience. 
We begin by recalling how the potential is written in terms of diagonal and off diagonal 
elements (eq. (^)). 

U(d, Q) = C int = fiu 2 Tr [-2(4 - d{){d ai - d aj )Q Uj Q Uj + 2« - di)(d M - d b ^)Qu 3 Q a31 

+ 4(4 - d{)Q Uj [Q a , Q b f + [Q a , Off] (48) 
The classical equations of motion are then, 

dl = u 2 



-8(4 - di)Q bij Q" ij + 8(d bl - d bj )Q Uj Q aji + 8Q bij [Q a , Q b f + i[Q b , [Q a , Q 

(49) 

Qtj = ^ - di) 2 - (4 - di)(d u - d b i))Q% + 4(4 - d{)Ql 3 - 8(4 - di)[Q a , Q b } tJ 

(50 

Now we will consider how each matrix element moves in an effective potential created by the 
averaged motions of the other elements. To see this we make a mean field approximation, 
good at large N. We assume that the statistical averages satisfy relations consistent with 
the symmetry of the theory. This gives us 

< Q^Qf >= q 2 5 ab (5 ik 5ji + 8 u 5 jk ) (51) 

< (4 - 4)(4 - 4) >= r 2 5 ab (5 lk ^ 1 - 5 a V k ) (52) 

< Qaij >=< 4 >=< d Q >=< Q 3 >=< d3 >= (53) 

We assume also that averages of four matrix elements factor into pairs of averages of two in 
all ways. Our goal will be to solve for the value of q when the off diagonal elements are in 
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a thermal bath. We now write out the effective potential for the matrix elements moving in 
the averaged fields of the other elements, to quadratic order. This gives us, 

<U>=^Q mi Q a - + ^K-dtf (54) 

with 

n 2 Q = A(d-l)uj 2 [(N -l)q 2 + 2r 2 } (55) 

n 2 d = A(d-l)uo 2 q 2 (56) 

Thus, each off diagonal element moves in a harmonic potential created by the averaged 
values of the other elements. The diagonal elements are, to leading order a system of points, 
each connected to all the others by springs with the same spring constant. We will make 
no assumption about the statistical distribution of the diagonal elements. We will see that 
this is consistent so long as p is chosen so that the diagonal spring constant fid vanishes as 
N — ► oo. We will assume that the Qaij's are in thermal equilibrium with each other. This 
implies that at temperature T, 

7^0 = T 
2 Q 4 

which tells us that 



= T (57) 



w^- = [N - l)q,+2r2q2 (58) 

Now let us assume that we can neglect the term in r 2 q 2 for estimating the order of magnitudes 
of quantities and the dependence on N. This will be consistent so long as r 2 < Nq 2 . We 
then arrive at the estimate that for large N, 



T \ 1/4 , . 1 



q 



iV 1 / 4 \8(d — \)[iio 2 ) N 1 ^ 



* 1/4 -i+r- (59) 



We see that for p — 1, it is true that Nq 2 is of order unity, so this is consistent also with r 
being of order unity as N — > oo and T — > 0. Note that we now have 

Q 2 Q = Au 2 (d - 1) [2r 2 + ViN 1 ^} (60) 

We see that for Qq not to diverge as iV — > oo we must have p > 1. We then see that p — 1 
corresponds to a critical point at which Qq is fixed as iV —>■ oo. We also see from (|59|) that 
this corresponds to q ~ 1/y/N. We see also that, 

O 2 n 2 1 

Vl 2 Q Nq 2 + 2r 2 N + 2 J 1 J 

Thus we see that if we choose p = 1 so Qq is fixed in the limit iV — > oo then in the same 
limit Qd vanishes, so the diagonal elements remain free in the mean field approximation. 
From now on, unless otherwise specified, we will take p = 1. Now we would like to estimate 
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the effects on the matrix elements of the fluctuations around the thermal averages for small 
temperature and large N. Let us consider for example the contribution of the last term in 
eq. p9| ) to the motion of d l a . We have, neglecting the other terms, 

tf a = ±u 2 Y\QbAQ a ,Q h ]h (62) 

bj 

Counting the commutators, this is a sum of dN 3 terms, for each a and i. These add with 
random signs. Each impulse of a fixed sign lasts an average time of fig 1 , as this is the time 
over which the values of the Qaij's oscillate. The result is that d l a has a random Brownian 
motion on top of the free motion given by the N — > oo limit. To estimate the diffusion 
constant that results note that if, in a time At the average displacement resulting from the 
random forces is Ad, the diffusion constant is 

"< = ~aT (63) 

If ^totai is the total acceleration given by the sum of the random forces over the time At we 
have Ad = \a to tai{.At) 2 . Thus, we have, 

v d =~a 2 total (At) 3 (64) 

For a tota i we may take 

atotai = ^ 2 g 3 v / iV^ (65) 

because we are adding N 3 d terms with random signs and average magnitude uj 2 q 3 . Taking 
At = Qq 1 we have for large N, 

ud q 3 N 3 ' 2 udt 3 l A 1 

Vd ~ 4(d - If' 2 [1 + |J]3/2 - 4(d _ l)3/ 2iV i(p-l) [1 + ((>(,) 

So we see that if we pick p — 1, and so long as 2r 2 is of order one, the diffusion constant 
for the diagonal elements goes to a limit which is iV-independent and hence T independent 
as iV — ► oo and T — > 0. Under these same conditions we can show that as iV — ► oo these 
terms make the dominant contribution to the random forces on the diagonal elements coming 
from eq. (f4Tj|). A similar analysis starting from eq. ( |5T^ ) allows us to estimate the diffusion 
constant for the Q's coming from the random forces to be, 

r 2 t 

The result is that for p > 1, vq vanishes in the limit N — > oo. Now a classic result of 
random matrix theory is that if the matrix elements of a matrix undergo Brownian motion, 
so do its eigenvalues. In a case like ours it is clear, as the eigenvalues will be close to the 
diagonal values, as the off diagonal values vanish as iV — >• oo. But we have to be careful 
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about the contributions from higher order terms in perturbation theory. By making use of 
the standard formula, 

A? = *? + E#% + - (68) 



d a — d a 



we can show that the diffusion constant for the eigenvalues is given by 



2 



v \ — v d^ — 2 h ... w w 



+ 



4(^_ 1)3/2^1(^-1) N9p/4+3/4 



(69) 



This tends to v& as N — > oo. Thus, we see that the dominant contribution to the Brownian 
motion of the eigenvalues comes from the Brownian motion of the diagonal elements. These 
in turn are fluctuating because they are perturbed by their interactions with the off diagonal 
elements, which are moving in a harmonic potential, created by their averaged values at 
finite temperature. The result is that a randomness is introduced into the motions of the 
eigenvalues, coming from the interactions of the diagonal elements with a very large number 
of random variables, which are the off diagonal elements. This then illustrates that idea 
that a local degree of freedom can have its motion randomized by interaction with a large 
number of non-local degrees of freedom. In the next section we will see that this may result 
in behavior that for large N is indistinguishable from that predicted by the Schroedinger 
equation. 



6 Derivation of the Schroedinger equation 

We are now ready to derive the Schroedinger equation for the eigenvalues of the matrices. 
As we described above this is a three step process, 

• STEP 1 Formulate the statistical variational principle for the matrix model. 

• STEP 2 Make assumptions about the statistical ensemble. In particular, we assume 
that the model is in an S'-ensemble, heat it to finite temperature T and then study the 
large N limit with Tml/N. 

• STEP 3 Derive an effective statistical variational principle for the eigenvalues by 
averaging over the variational principle of the matrix elements and show that when 
N — > oo this is equivalent to Schroedinger quantum theory for the eigenvalues. 
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6.1 STEP 1: The statistical variational principle for the matrix 
model 

We begin by denning an ^-ensemble for the matrix elements. That is, we begin with the 
variational principle 



I[p,S] = Jdt f(dd)(dQ)p(d,Q,t) 



(70)" 

where U(d, Q) is the interaction term C int given by (|^). 
6.2 STEP 2: Physical assumptions 

We than will state the physical assumptions we make concerning p and S. These are assumed 
only to hold to leading order in 1/N 

• The Q system is in a distribution that is to leading order in 1/N statistically inde- 
pendent of the distribution of the eigenvalues. This means that to leading order the 
probability density factorizes 

p(d,Q)=p d (d)p Q (Q) + 0(l/N) (71) 

• The Q subsystem is in thermal equilibrium at a temperature T. So we have 

PQ(Q) = \e-^' T (72) 



where H(Q) is the hamiltonian corresponding to the Q system alone 



H(Q)=p 

and 



£(^) 2 -u> 2 [Q a ,Q*][Q fl ,Q bl 

aij 



(73) 



Z = j dQe- H ^' T (74) 
As a result of these assumptions our variational principle reads, 

I[p d ,S,T\ = ! dt f(dd)(dQ)p d (d)p Q (Q) 



(75) 
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6.3 STEP 3: Derive an effective variational principle for the eigen- 
values 



Now we want to derive an effective variational principle to describe the evolution of the 
probability distribution for the eigenvalues. We will do this by averaging the variational 
principle (|75|) over the values of the matrix elements, and then extracting the leading behavior 
for large N and small T. We begin by inserting the factor unity in the form 



QiiQi; 



j U 3 



+ 



(76) 



Thus, we have, 



I\p d ,S,T} 



1 

Z 



dt / dddQ / d\5 A? - df - 



QijQji 

K - A? 



+ 



pde 



-H{Q)/T 



S + 



1_ 5S_, 
2p [ 5df 



+ + 



1 5S 



+ U(d,Q) 



(77) 



Now we would like to integrate over the gPs, which will express the action in terms of only 
the A's and Q's. However, before doing this we need to take into account that while the to 
the order we are working, the gPs and the A's will be undergoing Brownian motion because 
the diagonal elements are moving in a random potential given by the values of the Q's. 
There are also additional contributions to the diffusion constant coming from the terms in 
Q that contribute to the eigenvalues at higher order . So we will have to be careful about 
the definitions of the velocities. In particular, we will have to recall that in the theory of 
stochastic processes the limits which define time derivatives are taken after the averages over 
probability distributions, not before. So as before we must write 



J dt j dddQp d {d)p Q {Q)^-{^) 2 = J dtdddQ Pd {d)p Q {Q))p{V{d) 

'd?(t + At)-dm 2 \ 



2/i 

dt hm^J dddQpd(d)pQ(Q)p 

dt A™oJ ddd Qp^ d )pQ^Q)^ 



At 2 J 

f d?(t+At)-dmr 

At 2 



(7E 



' d«{t) - d«{t - At)f 
At 2 



Note that the last equation follows trivially, for smooth motion, but it will have non trivial 
consequences once we have averaged over the Q's because the result for large N is to induce 
Brownian motion for the off diagonal elements and eigenvalues. Now we perform the integral 
over the gPs. It is useful to write 

dl = A? + AA? (79) 

where 

QijQji 



AA«(Q,A) = -E 



A? 



+ 



(80) 
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has to be treated as a stochastic variable, taking into account its dependence on the Q's 
which are themselves fluctuating due to the assumption that they are in equilibrium in a 
potential. We then have, to leading order in 1/N, 



I[p d ,S d ,T] = JdtJ dXp d (X,t) J dQ PQ (Q) 

+ K.E. (81) 



1 .5S(X + AA,Q). 2 



S(X + AA, Q) + — ( v — ^' Y + U(X, Q) 
where to leading order, the kinetic energy terms for the d's have become, 

n \ f(\t(t + At) - xmr\ + ({m) - xt{t - At)r 



K.E. = J dtlimj dXp d (X,t) J dQp Q (Q)^ 



At 2 J \ At 2 

(82) 

We now are ready to integrate over the Q's. The key point is that the dependence of the 
A's on the Q's through a sum of a large number of independent terms, J2j v 7 ,^' , as well as 

i 3 

the coupling of the A's with the Q's coming from the terms in U(X, Q) turns the A's into 
stochastic variables, described by a stochastic differential equation of the form 

DX a i = 6 ? °(A, t)dt + AA° At > (83) 

DX a i = b* a {X, t)dt + A*Xi At<0 (84) 

with 

< AA^AAj >= S ab 5 il u x dtdt > (85) 

< A*A l °A*A$ >= -5 ab 5 i:j u x dt dt<0 (86) 

< F(X, Q) >= J dQ PQ (Q)F(X, Q) (87) 

We note that we can use the value of v\, ( |69| ) computed in the last section as all the 
assumptions we made there have been carried over here, so long as we work to leading order 
in 1/N with T scaled as T m 1/N. We also have, from the Focker-Planck equations that the 
current velocity is 

v?W = \w + br) (88) 

while the osmotic velocity is 

<(A) = ^?-C) = ^^ (89) 



Here the brackets mean 



^From these we can derive, 

hm/ dQp d (X, t)p Q (Q)~ ( {K{t + m)2 ) = p d (X,tW(X 7 tf + NC] (90) 
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and 



hmj dQp d (X,t)p Q (Q)- UU ^ 2 U- )=p d (\,t)[b?(\,t) 2 + NC] (91) 



where C is the infinite constant we defined in eq. (|33|) . 

To go further we need to define the effective Hamilton- Jacobi function for the eigenvalues. 
We define 

S x (\) = J(dQ)p Q (Q)S(X,Q) (92) 
We now show that, to leading order in 1/N, 

= (93) 

Consider the probability conservation law that follows from the statistical variational prin- 
ciple that defines the dynamics of our matrix model, eq. (f7C|). 



p(d,Q) = -- 
p 



s Mi<Q) s J^m )+ j_ mQ) s jm}. ) 



(94) 



_3d a i $d a . t oQaij $Qaij 

But using (|7T| ) and (|72|) we have that 

p(d,Q)= Pd (d)p Q (Q). (95) 

We also have, by the same assumptions, since a thermal distribution is stationary and has 
no current velocity, 

v a HQ) = - 6 -^^ = 0(l/N) (96) 
Thus, we have, integrating over the Q's, 

Pd(d) = --JrMQJr I dQ PQ (Q)S(d,Q)) + 0(l/N). (97) 



p 5d ai 5d, 

To leading order we can replace everywhere the dependence on d a i with dependence on A 
since the terms by which they differ are also higher order in 1/N. Thus we have 



P*W = --7T-(^( A )^T^) + °(VA0 (98) 
p b\ ai 0A a 



But by (UK) we must have 

OAai 

This establishes eq. (|93|). With this result we have the key relation that, 

+ n = + «>) = M)° + ^ paep) > . (ioo) 
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We also define 



and 



E Q = J dQ PQ {Q) 



l t 8S Q {Q), 2 + f rf TrUQ? 
2 



f nSEW - A D 2 = / ^p q (Q)^(a,q) 



(101) 



(102) 



aij 



We can estimate that Eq = TN(N - l)/4 « Npu 2 so this is a divergent constant in the 
limit. The result is 



I[p d ,S,T] = fdtj d\[p d (\,t)S x -H e "(S x ,p d ,T)} 
where, the effective hamiltonian for the eigenvalues is 



(103) 



H*"(S x ,p x ,T)=p d (\) 



E' Q +< 



&8S X (\)\ p V \ [8\np x (\y 



5\? 



+ 



5\? 



^E(a?-a?) 2 

aij 



(104) 

where E' Q = Eq + NpC contains both infinite constants. The resulting equations of motion 
are 



E ' n ,S' A - J- f^^l + ^E( A ? " A ^ )2 + """ = 



2^V $K 
and the current conservation equation 



at] 



px = —d ai p x (d ai S x ) 
p 



(105) 



(106) 



The so-called "quantum potential" is given by 

quantum = ^2 f (^^1 ^ + l^(p A ^ In 



;• pi4^=^'\^) (107) 

VPa(A) 



These we recognize as the real and imaginary parts of the Schroedinger equation, when we 
write 

^(X,t) = ^p- X e s ^ h (108) 

with 

dt 3 / 2 

h = pv x = /iw— — - vq/9 (109) 



4(rf- l) 3 /2 



So, finally, we have in the limit N — > oo, 



«7i 



ft 2 <5 2 



2/i5(A 



A^d ^Ad \a\2 i en/ 

Q 



*(A,f) 



(110) 
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Finally, by repeating the argument of section 3.4 we can show that the conserved energy of 
the original theory splits into two pieces, 



H = H* + E' n 



where 



ti 2 5 



2A* *(A?) 



a\2 



IP. 



(Ill) 



(112) 



aij 



Since E'q is an infinite constant the result is that if*, which is the quantum mechanical 
energy, is conserved as iV — > oo. We can then renormalize the wavefunctional so that 



* r (A) = e lE Q t/h V(\) 



(113) 



Finally, we note that as Q 2 , ~ 1/N the eigenvalues become free in the limit iV — > oo. Thus, 
when iV — > oo the probabilities evolve according to the free Schroedinger equation, 



th 



d* P (A,t) 
dt 



h 2 5 2 



(114) 



7 Conclusions 



What we have shown here may be summarized by saying that matrix theory may not only 
give rise to string theory, and hence gravity, it may also give rise to quantum theory, in 
the sense that the quantum evolution of the eigenvalues may appear, at large N, to be a 
consequence of the classical statistical physics of the matrices. This means that we may be 
able to solve the daunting conceptual problems of quantum theory by means of a simple 
physical hypothesis: that the theory of gravity and hence spacetime arises from a non-local 
background independent theory in which geometry initially plays no role and the physical 
degrees of freedom represent relational rather than intrinsic properties. There remain many 
open questions. A short list is: 

• Is it possible to extend these results to the actual supersymmetric models that arise in 
string theory^? 

• Is it possible that the existence of dualities that connect certain quantum field theory 
observables to the classical limit of string theory are related to the fact that the classical 
matrix theory can in a certain limit reproduce a quantum theory? 

• How does Lorentz invariance and relativistic causality arise in the matrix models which 
describe relativistic string and membrane theories, and how is this compatible with the 
non-local dynamics of the off diagonal elements? 

4 This is in progress with Stephon Alexander JTo| . 
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Can this be extended to truly background independent matrix models, such as the 
cubic matrix models [12] and the matrix models which have been developed for spin 
I? 



foams [1X3 



Are there any practical experimental predictions that follow from these theories? 
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